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1  introduction 


Background 

The  Defense  Nuclear  Agency  funded  the  Explosion  Effects  Division  (EED) 
of  the  USAE  Waterways  Experiment  Station  (WES)  to  develop  a  large-bore 
(4-ft-diameter)  vertical  gas  gun.  The  primary  purpose  of  the  gas  gun  is  to 
simulate  ground  shock  environments  in  a  manner  suitable  for  testing  ground 
shock  measurement  devices.  An  additional  purpose  for  the  gun  is  to  obtain 
dynamic  material  properties  data  from  relatively  large  geologic  samples. 


Description  of  the  DN A/WES  Ground  Motion  Test 
Facility 

Design 

The  facility  was  designed  over  a  24-month  period  beginning  in  1988  (see 
White  et  al.  1991).  An  artist’s  rendering  of  the  facility  is  shown  in  Figure  1, 
and  a  schematic  of  the  gun  is  shown  in  Figure  2.  The  gun  consists  of  a  large 
annular  pressure  vessel  surrounding  a  vertical  barrel.  A  series  of  orifices  are 
machined  in  the  upper  barrel  wall.  These  allow  the  compressed  air  from  the 
vessel  to  expand  into  the  barrel.  When  the  gun  is  in  the  cocked  position,  a 
projectile,  with  o-rings  placed  at  the  top  and  bottom,  straddles  the  orifices, 
and  prevents  the  compressed  air  from  being  released  into  the  barrel.  The 
projeaile  is  held  in  place  by  a  quick-release  trigger  mechanism.  A  water 
reaction  mass  fills  the  top  portion  of  the  barrel  above  the  trigger  mechanism. 
This  mass  is  used  to  reduce  the  reaction  of  the  gun  to  the  pressure  drivmg  the 
projectile.  The  bottom  of  the  barrel  may  be  sealed  with  a  di^hragm  to  allow 
a  partial  vacuum  to  be  created  in  the  barrel  seaion  below  the  projectile. 

(Note;  none  of  the  tests  rqxirted  here  utilized  a  vacuum  in  the  lower  portion 
of  the  barrel.) 

Operation 

To  frre  the  gun,  the  projectile  is  released.  The  weight  of  the  projectile 
causes  it  to  move  downward.  As  the  top  o-ring  clears  the  orifices,  the 
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Figure  2.  Schematic  of  the  4-ft-diameter  gas  gun 
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compressed  air  expands  into  the  barrel.  The  incoming  air  simultaneously 
drives  the  projectile  downward  and  the  water  reaction  mass  upward.  Theoret¬ 
ically,  the  mass  of  the  water  can  be  adjusted  so  that  the  "bottom”  of  the  mass 
of  water  will  exit  the  top  of  the  barrel  at  the  same  time  that  the  projectile 
clears  the  bottom  of  the  barrel.  This  would  maintain  maximum  pressure  on 
the  projectile  for  as  long  as  it  was  in  the  barrel.  The  maximum  operating 
range  is  2,070  kPa,  which  produces  a  projectile  velocity  of  71  mps. 


Objectives 

The  primary  objectives  for  all  tests  were  to  develop  safe  procedures  for 
testing  and  to  determine  the  operating  characteristics  of  the  gun  (e.g.,  determi¬ 
nation  of  projectile  velocity,  far-field  ground  motion  and  nuisance  noise  levels 
as  a  function  of  vessel  pressure).  The  performance  aspects  of  testing  with  the 
4-ft  gun  with  regard  to  the  primary  objectives  have  been  reported  by  White 
(1993). 

The  primary  objectives  of  the  development  program  for  the  gas  gun  could 
have  been  met  without  including  testbeds  on  each  test.  However,  to  take  full 
advantage  of  the  gun,  testbeds  were  included  when  possible.  This,  in  turn, 
led  to  a  secondary  set  of  objectives  being  associated  with  these  testbeds. 

These  objectives  included  the  initial  dynamic  testing  of  new  ground  shock 
instruments  for  use  in  the  HUSKY  JAGUAR  and  the  Seismic  Hard-Rock  In 
Situ  test  programs,  and  the  development  of  testing  procedures  and  equipment 
for  routine  use  widi  the  gun.  The  inclusion  of  a  target  on  a  test  never  became 
the  driving  force  of  the  gas  gun  program;  however,  its  role  increased  in  im¬ 
portance  as  testing  progressed.  Of  the  twenty-two  tests  conducted  under  this 
program,  twelve  of  these  tests  included  an  instrumented  testbed. 


Scope 

This  rqwrt  describes  the  testbed  construction  and  test  procedures.  Stress 
and  velocity  histories  are  presented  from  each  test.  Analysis  of  the  data 
includes  determination  of  the  propagation  velocity  of  the  loading  wave  for  the 
two  types  of  sand  used  as  a  backffll.  Two  methods  for  determining  the  propa¬ 
gation  velocity  of  the  rarefaction  wave  are  discussed.  Two  methods  for 
analyzing  impaa  angle  and  planarity  are  also  presented.  Conclusions  are 
drawn  on  the  results  of  the  development  program  to  date. 
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2  Experiment  Description 


Testbed  Construction 

Twelve  tests  were  conducted  using  a  sand  testbed  as  a  target  beneath  the 
gun  during  the  period  addressed  by  this  report  (January  1990  •  December 
1992).  Nine  of  the  tests  used  masonry  sand  and  three  of  the  tests  used 
Socorro  Plaster  sand  as  the  backfill  material.  The  disposable  target  container 
was  a  1.37-m  diameter  by  0.76-m  high  corrugated  12-gauge  steel  pipe  with  a 
wdded  12-gauge  steel  bottom. 

The  target  for  each  test  was  constructed  by  placing  four  lifts  (typically 
lS-20  cm)  of  sand  in  the  steel  container.  Each  lift  was  compacted  using  a 
vibratory  compactor  to  obtain  a  maximum  density.  The  wet  density  and 
moisture  content  were  measured,  using  a  nuclear  densitometer,  at  two 
locations  in  each  of  the  three  lower  lifts.  For  several  tests,  the  moisture 
content  was  also  measured  by  weighing  samples  of  the  sand  before  and  after 
drying  in  a  microwave  oven  (called  the  oven  dry  method).  The  last  test.  Test 
22,  was  constructed  to  meet  specifications  for  density  and  moisture  content. 
The  construction  methods  used  for  the  othw  tests  were  designed  primarily  to 
produce  a  uniform  testbed,  as  opposed  to  specific  densities  or  moisture 
contents. 

The  time  required  to  build  a  testbed,  once  all  materials  were  on  hand 
(gages,  sand,  vibrator,  microwave  oven,  tools,  etc.),  was  approximately  one 
day  per  test.  An  efficient  m^od  for  testing  with  the  gun  was  to  pr^are 
several  samples  at  a  time,  since  objectives  for  subsequent  tests  were  not 
contingent  on  results  from  previous  tests.  For  this  series,  two  or  three 
testbeds  were  prq)ared  (in  a  sheltered  area)  on  successive  days.  Using  a  fork 
lift,  they  were  moved  to  the  facility  as  required  for  testing.  Approximately 
two  to  three  days  were  required  for  conducting  a  test. 


Typical  Instrumentation  Layout 

The  total  d^th  of  each  target  was  61  cm.  Two  depths  in  a  target  were 
typically  instrumented.  These  two  depths  (15.2  cm  and  30.5  cm)  included  a 
variety  of  the  stress  and  velocity  instruments  (Tables  1  and  2).  The  conunon 
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Table  2  I 

Velocity  Gages  Used  for  Tests  in  Sand  | 

1  Inatniinant/Caniatar 

Common  Name  1 

1  Log  Shook-laolated  Accelerometer  canieter 

AV-L  1 

1  Wadge  Shock-laolated  Accalarometar  carxeter 

AV-W  1 

1  Micro  Hard-Mounted  Accelerometer  canieter 

AV-M  1 

1  HIFI  Accelerometer  cartieter 

AV-N  or  SHMA  1 

1  Street  and  Velocity  gage 

S&V  or  CSHMA  1 

1  Sunburst  Recovery,  Inc.  Constant  Pitch  Electromagnetic  gage 

SC 

1  Sunburst  Recovery,  Irtc.  Variable  Pitch  Electromagnetic  page 

SV 

name  used  to  describe  each  instrument  is  also  listed  in  the  tables.  A  rq)resen- 
tative  example  of  the  gage  layout  is  presented  in  Figure  3.  Gages  were  placed 
in  die  correct  lift  as  the  testb^  was  being  built.  The  instruments  located  at 
the  30.5  cm  dqith  were  typically  placed  at  a  radius  of  15.2  cm,  and  the  in¬ 
struments  at  the  15.2  cm  d^th  wm-e  placed  at  a  radius  of  30.5  cm  (see  Fig¬ 
ure  4).  For  some  tests  a  gage  was  placed  in  the  cmiter  of  the  testbed  at  one  of 
the  two  levels.  A  short  length  of  cable  was  routed  from  each  gage,  through 
the  testbed,  and  terminated  just  outside  of  the  target  container.  The  gages 
used  in  eadi  testbed  and  their  locations  are  listed  in  Appendix  B. 

Diagnostic  measurements  were  included  for  each  test.  Two  pressure  gages 
were  located  in  the  reservoir  of  the  gas  gun.  One  served  as  a  direct  readout 
gage  used  to  monitor  system  pressure,  prior  to  firing,  at  the  control  panel  for 
the  gun.  The  second  gage  was  used  to  monitor  the  pressure  drop  in  the  vessel 
during  a  test.  An  accelerometer  was  mounted  on  the  base  of  the  gun  to 
monitor  acceleration  levds  in  the  vertical  direction.  Six  piezoelectric  pins  at 
the  bottom  of  the  barrel  of  the  gun  were  used  to  measure  the  projectile 
velocity  and  planarity  as  the  forward  face  exited  the  barrel.  A  triaxial  seismic 
station  and  microbarogrt^h  gage  were  fielded  at  three  different  ranges  on  Sev¬ 
a’s!  tests  to  monitor  far-freld  ground  motions  and  airblast.  Sound  pressure 
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DEPTH  =  15.2  CM 
RADIUS  =  30.5  CM 


0°  HRSE-1 


Figure  3.  Typical  layout  for  testbed  instrumentation 
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Figure  4.  Cross  section  of  the  gas  gun  projectile  and  container  of  target  material 
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level  measurements  were  made  at  five  ranges  on  most  tests,  and  a  video 
camera  was  used  to  photograph  each  test.  The  results  of  die  diagnostic  mea¬ 
surements  have  been  rqxirted  by  White  (1993). 


Testing  Procedure 

For  these  tests,  die  projectile  was  loaded  into  the  gun  prior  to  moving  the 
tesdied  beneath  the  barrel  of  the  gun.  The  projectile  was  raised  into  its 
cocked  position  via  a  chain  hoist  attached  to  a  tripod  at  the  top  of  die  barrel. 
Two  chains  were  used  to  lift  the  projectile,  and  two  were  used  as  saf^  chains 
during  the  loading  process.  The  two  safety  chains  remained  in  place  after  the 
projectile  was  locked  into  position  with  the  quick-release  trigger  mechanism. 

The  configuration  of  the  projectile  varied  d^iending  on  the  vessel  pressure 
planned  for  the  test  (see  Figure  4).  As  test  levels  increased  (higher  vessel 
pressures  and  projectile  velocities),  thicker  pieces  of  foam  were  required  at  the 
base  of  the  projectile  to  absorb  the  impact  energy,  and  thus  protect  the  car¬ 
riage  from  damage.  The  mass  of  the  projectile  varied  slighdy  as  a  result  of 
changing  the  foam  thickness.  The  mass  of  the  projectile  for  most  tests  was 
1,482  kg.  The  vessel  pressure  for  these  twelve  tests  varied  from  34S  to 
1,379  kPa  and  the  projectile  velocity  (as  the  leading  edge  exited  the  barrel) 
varied  from  27.9  to  S9  mps.  The  location  of  the  target  was  typically  120  cm 
below  die  end  of  the  barrel.  The  test  level,  projectile  mass,  and  projectile 
velocity  for  each  test  is  listed  in  Table  3. 

Sand  was  placed  over  the  earthen  portion  of  the  trench  bottom  located 
direcdy  beneath  the  barrel  of  die  gun.  Several  timbers  (4-in.  by  6-m.)  were 
placed  in  die  sand  base  to  support  the  sanqile.  The  sample  was  moved  into 
position  beneath  the  barrel  of  the  gun  via  a  trolley  that  ran  along  channels  cast 
into  the  concrete  walls  on  either  side  of  the  trench.  The  sample  was  centered 
beneadi  die  gun  using  a  plumb  bob  suspended  from  crosshairs  at  the  bottom  of 
the  barrel.  Final  leveling  of  the  testbed  was  achieved  using  a  carpenter  level 
to  guide  placement  of  shims  betwem  the  sted  bottom  of  the  sample  and  the 
timbers.  An  advantage  of  supporting  the  sample  with  timbers  is  that  it  essen¬ 
tially  provides  a  "free"  surface  at  its  base.  The  shock  wave  reflecting  off  this 
free  surface  may  be  seen  in  the  stress  wave  forms,  and  thus  may  be  used  to 
investigate  the  unloading  wave  speed  of  the  testbed  material. 

After  the  testbed  was  positioned  beneath  the  gun,  instrument  cables  were 
spliced  to  cables  running  to  the  Control  Trailer.  After  all  electrical  checkouts 
were  completed,  die  safety  chains  on  the  projectile  were  removed.  A 
1-mm-thick  fiberglass  diaphragm  and  plastic  liner,  used  for  containing  the 
water  reaction  mass,  were  placed  above  the  quick-release  trigger  mechanism 
in  the  upper  portion  of  the  barrel.  The  depth  of  the  water  reaction  mass  for 
all  tests  was  1.3S  m,  which  corresponds  to  a  mass  of  1,570  kg.  After  placing 
the  reaction  mass,  the  firing  system  for  the  trigger  mechanism  was  enabled  by 
supplying  pressurized  gas  to  the  system  and  by  connecting  battery  power  to  a 
solenoid  valve  used  to  control  the  flow  of  gas. 
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Table  3 

Teat  Pressure  Level.  Projectile  Mess,  and  Projectile  Velocity  for 
Tests  in  Sand 


Teat  No. 

Veaaal  Praaaure 

kPa 

Projectile  Maaa 

ks 

Projeetila  Velocity  1 

"4-  - 1 

6 

34S 

1,4S2 

27.9 

7 

617 

1,482 

37.1 

8 

S17 

1,482 

35.5 

10 

862 

1,489 

47.2 

11 

1,014 

1,475 

50.3 

12 

662 

1,482 

41.8 

13 

1,027 

1,482 

52.8 

14 

1,207 

1,482 

55.7 

IS 

1,379 

1,482 

59.0  1 

16 

841 

1,482 

46.3'  1 

17 

1,1 6S 

1,482 

54.9'  1 

22 

862 

1,472 

46.8  1 

- - - 

[  *  Calculated  v«|y  (fw  expariment  data) 


Instrument  recording  and  remote  operation  of  the  gun  were  performed  at 
the  Control  Trailer,  located  q)proximately  SO  m  away.  A  picture  of  the  re¬ 
cording  s^p  is  shown  in  Figure  S.  The  control  panel  used  for  operating  the 
"till  systmn"  for  the  gun  is  shown  in  Figure  6.  llie  fill  syston  (illustrated  in 
Figure  7)  is  comprised  of  three  components;  the  pressure  system,  the  vacuum 
system,  and  the  firing  system.  None  of  the  tests  rq)orted  here  used  a  vacuum 
in  foe  lower  section  of  foe  barrel.  The  pressure  vessd  was  filled  by  a  hi^- 
pressure,  high-volume  air  compressor.  Alter  foe  desired  testing  level  was 
reached,  foe  mr  compressor  was  stopped.  The  time  required  to  pressurize  foe 
gun  for  these  tests  varied  between  13  and  54  minutes.  A  master  switch  was 
then  used  to  ensure  all  components  of  foe  pressure  system  were  in  their  proper 
state;  i.e.,  all  valves  closed  and  foe  air  compressor  shut  off.  The  master 
switch  also  provided  power  to  foe  firing  system.  After  sounding  a  warning 
sirm  and  a  final  check  with  instrumentation  personnel,  foe  gun  was  fired  by 
energizing  foe  4-way  solenoid  valve  in  foe  firing  system.  This  allowed  gas 
pressure  to  activate  foe  quick-rdease  trigger  mechanism,  which  released  foe 
projectile.  On  average,  approximately  450  msec  were  required  for  foe  triggo- 
medianism  to  release  foe  projectile.  The  additional  time  required  for  foe  pro¬ 
jectile  to  travel  foe  length  of  foe  barrel  was  anywhere  between  1(X)  and 
400  msec,  depending  on  foe  test  level  and  len^  of  foe  projectile. 
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3  Experiment  Data 


Testbed  Material  Data 

The  wet  density  and  moisture  content  of  the  test  materials  were  measured 
at  two  locations  in  each  of  the  three  lower  lifts  of  a  testbed  using  a  nuclear 
densitometer.  The  accuracy  of  the  nuclear  densitometer  for  measuring 
moisture  content  is  questionable,  however,  because  the  source  and  receiver  are 
both  in  the  base  of  the  gage,  and  experience  has  shown  that  the  moisture 
content  is  only  measured  to  a  d^th  of  approximately  S  cm.  For  that  reason, 
the  moisture  content  was  also  measured  for  several  tests  by  weighing  samples 
of  the  sand  before  and  after  drying  in  a  microwave  oven  (oven  dry  method). 

A  summary  of  the  average  wet  density  and  moisture  content  of  the  entire 
testbed  is  presented  in  Table  4.  Presented  in  Appendix  C  are  profiles  of  the 
wet  density  and  moisture  content  for  each  tesfoed.  The  values  presented  in 
the  appendix  r^resent  the  average  of  two  measurements  at  each  d^th. 

Who'e  moisture  content  measurements  were  made  with  both  the  nuclear  densi¬ 
tometer  and  die  ovoi  dry  method,  the  oven  dry  m^od  (as  the  more  accurate) 
is  presented.  Considering  the  average  moisture  content  of  the  entire  testbed, 
the  oven  drying  method  indicated  values  anywhere  between  0.4  and 
1.8  percent  higher  dian  those  determined  using  the  nuclear  densitometo^. 


Active  Measurement  Data 


Stress  and  impulse  measurements  from  each  test  are  presented  in  Fig¬ 
ures  D1-D12,  contained  in  Appendix  D.  Acceleration  records  and  the  inte¬ 
grated  velocity  and  displacement  records  are  presented  in  Figures  D-13 
through  D-24.  Data  return  from  the  tests  was  very  good.  The  average  peak 
stresses  produced  in  the  test  materials  varied  between  7  and  27  MPa.  The 
average  peak  particle  velocities  varied  from  18  to  34  mps. 
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Table  4 

Average  Wet  Density  and  Moisture  Content  for  Sand  Testbeds 

Teat  Me. 

Sand  Type 

Aeaa^  1MM  DanaMy 

i 

1 

il 

6 

Meaonry 

1.61 

4.9<" 

7 

Meeorwy 

1.59 

3.6“ 

8 

Maaorwy 

1.61 

5.2“ 

10 

Meaofwv 

1.60 

3.5“ 

11 

Meaonry 

1.62 

3.8“ 

12 

Maaonry 

1.58 

3.0“ 

13 

Maaoitry 

1.58 

1.4“ 

14 

Maaonry 

1.59 

to 

i 

15 

Maaonry 

1.60 

3.1“ 

16 

Socorro  Platter 

1.69 

2.4“ 

17 

Socorro  Plaater 

1.70 

2.5“ 

22 

Socorro  Plaatar 

1.67 

4.0“ 

°°  Oven  dry  method 

Nuclear  dettaitometer 

Time>of-Arrival  Data 

The  propagation  velocity  of  the  shock  wave  in  sand  was  determined  by 
noting  the  time-of-arrival  (TOA)  at  known  dq>ths  widiin  the  testbed.  For 
eleven  of  the  tests  (all  but  Test  22),  instruments  were  located  at  only  two 
dq)ths  within  the  testbed  (15.2  cm  and  30.5  cm).  The  cento-line  of  each 
gage/canister  was  placed  at  the  designated  dq>th  (see  Figure  8).  After  the 
shock  wave  reach^  die  upper  edge  of  a  canister,  the  time  required  for  it  to 
traverse  to  the  sensing  element  was  n^ligible  compared  to  the  time  required 
for  the  shock  wave  in  sand  to  traverse  a  like  distance.  As  a  result,  the  rda- 
tive  "distance"  between  the  top  surface  of  die  various  instruments  was  used  to 
determine  the  "dqith"  over  which  the  shock  wave  propagated  between  the 
TOA  at  each  location. 

Each  test  in  sand  (excqit  Test  22)  included  a  Log  shock-isolated  accelw- 
ometer  (SIA)  canister  (AV-L)  at  the  15.2  cm  depth.  This  canister  had  the 
largest  thickness  of  those  used  in  the  tests,  and  for  planar  impacts  should  indi¬ 
cate  the  first  TOA.  Presented  in  Figures  9  throu^  19  is  the  TOA  at  each 
instrument  location  for  the  first  eleven  tests  in  sand.  The  x-axis  for  these  fig¬ 
ures  is  labeled  as  "distance,"  since  the  values  presented  are  die  relative  dis¬ 
tance  between  the  top  surface  of  a  Log  canister  and  other  gages/canisters  in 
die  testbed.  For  Test  22,  TOA  crystals  were  placed  at  the  surface  and 
diroughout  the  dqith  of  the  testbed.  The  TOA  data  presented  in  Figure  20  is 
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Figure  8.  Typical  target  geometry  for  tests  in  sand 


TOA  (MICRO-SEC) 
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Figure  9.  Time-of-arrival  vs.  distance  for  Test  6 
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Figure  10.  Time-of-arrival  vs.  distance  for  Test  7 
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Figure  1 1 .  Time-of-arrival  vs.  distance  for  Test  8 


TOA  (MICRO-SEC) 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

00 

<q 

CO 

cvT 

1— 

T~ 

1— 

1^ 

1^ 

1^ 

Chapter  3  Experiment  Data 


19 


Figure  12.  Time-of-arrival  vs.  distance  for  Test  10 
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Figure  13.  Time-of -arrival  vs.  distance  for  Test  11 
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Figure  14.  Time-of-arrival  vs.  distance  for  Test  12 
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distance  for  Test  13 


(MICRO-SEC) 


Figure  16.  Time-of-arrival  vs.  distance  for  Test  14 
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Figure  17.  Time-of-arrival  vs.  distance  for  Test  15 


MICRO-SEC 


Figure  18.  Time-of-arrival  vs.  distance  for  Test  16 
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Figure  20.  Time-of-arrival  vs.  distance  for  Test  22 


plotted  against  the  dq)tfa  (from  die  surface  of  the  testbed)  to  the  top  surface  of 
each  gage.  The  invme  slope  (P-wave  velocity)  of  die  best  linear  frt  dirough 
the  data  is  shown  in  the  figures  and  is  tabulated  in  Table  5  below. 


Tables 

Propagation  Velocity  of  Loading  Shock  Wave  for  Tests  in  Sand 


TaatNo. 

•andTypa 

Vaaaal  Praaaura 
kPa 

PreiactBa  Velocity 
mpa 

Wava  Propagation 
Volocity  mpo 

6 

Maaonry 

345 

27.9 

399 

7 

Maaonry 

517 

37.1 

491 

8 

Maaonry 

517 

35.5 

503 

10 

Maaonry 

862 

47.2 

328 

11 

Maaonry 

1.014 

50.3 

417 

12 

Maaoruy 

662 

41.8 

398 

13 

Maaonry 

1,027 

52.8 

337 

14 

Maaonry 

1,207 

55.7 

323 

IS 

Maaonry 

1,379 

59.0 

323 

16 

Socorro 

841 

46.3* 

280 

17 

Socorro 

1,165 

54.9’ 

309 

22 

Socorro 

862 

46.8 

'  Calculat«d  value  (no  expariniant  data) 


In  these  tests,  perfecdy  normal  impacts  are  unlikely.  Nonnormal  impacts 
can  result  in  apparent  inconsistencies  in  the  time-of-arrival  versus  dqith- 
to-top-of-canister  data.  For  instance,  in  most  tests  the  Wedge  SIA  canister 
was  placed  180**  from  the  Log  SIA  canister,  and  at  the  same  depth.  If  a 
nonnormal  impact  occurred  where  die  i^jectile  impacted  the  side  of  the  test¬ 
bed  near  die  Wedge  canister  first,  the  shock  wave  could  be  witnessed  by  the 
Wedge  canister  before  the  Log  canister.  This  is  seen  in  Test  1 1  (Figure  13). 
Notice  die  TOA  at  AV-W  occurs  prior  to  that  at  AV-L.  This  is  possibly  the 
explanation  for  several  cases  where  gages  at  slighdy  lower  depths  witnessed 
shock  arrival  prior  to  gages  at  higher  locations  in  the  testbed.  Other  explana¬ 
tions  for  this  might  include  accuracy  of  gage  placements  and  variability  in  the 
soil  conditions  across  the  diameter  of  the  testbed. 
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Determination  of  Propagation  Veiocity  of 
Rarefaction  Wave 


A  rarefaction  wave  was  generated  in  each  testbed  by  the  int^action  of  the 
shock  wave  with  the  free  surface  at  the  bottom  of  container  (and  the  sides  as 
well).  Two  methods  were  used  to  determine  the  propagation  velocity  of  the 
rarefrction  wave  (or  unloading  wave)  for  several  tests  in  sand.  Both  methods 
hinged  on  witnessing  the  TOA  of  the  unloading  wave  on  individual  stress 
records.  The  displacement  of  a  stress  gage  was  assumed  to  be  identical  to  the 
displacement  of  an  accelerometer  canister  located  at  the  same  d^th.  Discus¬ 
sion  of  these  results  is  limited  to  Test  22. 


Graphical  Method 

A  schematic  of  the  dqtths  of  the  stress  and  velocity  gages  fielded  in 
Test  22  in  their  original  position  and  at  a  later  time  is  presented  in  Figure  21. 
R^resentative  stress  records  from  the  IS.2  and  30.5-cm  d^ths  are  shown  in 
Figure  22.  The  TOA  of  the  rarefaction  wave  at  each  gage  location  is  noted  in 
the  figure.  The  displacement  of  the  accelerometer  canister,  at  each  dqith,  at 
the  arrival  time  of  the  relief  wave  is  measured  from  the  double-integrated  ac- 
cel^ation  history.  This  displacement  is  added  to  the  original  gage  position  to 
determine  die  location  of  the  gage  at  the  arrival  time  of  the  relief  wave.  The 
location  of  each  stress  gage  at  the  arrival  time  of  the  rarefaction  wave  is  plot¬ 
ted  in  Figure  23.  The  slope  of  the  linear  fit  through  this  data  indicates  an 
unloading  propagation  velocity  of  966  mps. 

If  the  loading  and  unloading  wave  velocities  were  constant  throughout  the 
dqith  of  the  testbed,  then  the  two  lines  r^resenting  the  best  fit  through  the 
loading  TOA  and  unloading  TOA  should  me^  at  the  bottom  of  the  testbed 
(l.e.,  at  the  61-cm  dqith).  The  nominal  stress  at  the  15.2-cm  d^th  as  record¬ 
ed  by  HRSE-1  was  27.1  MPa,  and  the  value  at  the  30.S-cm  dqiA,  as  record¬ 
ed  by  HRSE-2,  was  16.6  MPa.  This  difference  in  the  stress  levels  could 
cause  changes  in  the  loading  or  unloading  wave  velocities.  Other  potential 
reasons  for  the  non-intersection  of  the  lines  at  the  61-cm  depth  include  experi¬ 
ment  error  in  placing  the  gages  in  the  testbed,  or  uncertainty  in  the  absolute 
dq>th  of  the  testbed. 


Analytic  Method 

The  second  method  used  to  determine  the  propagation  velocity  of  the  rar¬ 
efaction  wave  in  the  soil  is  based  on  the  geometry  of  the  testbed  and  die  d^er- 
mined  loading  wave  velocity.  After  the  shock  wave  impacted  a  stress  gage,  it 
continued  until  it  impacted  the  free  surface  at  the  bottom  of  the  testbed.  The 
distance  over  which  the  shock  wave  traveled  at  the  loading  wave  speed  is 
known.  The  reflected  wave  then  traveled  back  though  the  testbed  at  the  un¬ 
loading  wave  speed.  The  distance  over  which  the  rarefaaion  wave  traveled 
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Figure  22.  Stress  measurements  at  the  two  instrumented  depths  from  Test  22,  indicating  TOA  of  the  loading  (T1  and  T2)  and 
unloading  {T3  and  T4)  waves 


(MICRO-SEC) 
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e  23.  Graphical  determination  of  the  unloading  wave  speei 


until  it  impacted  the  stress  gage  was  d^ermined  by  assuming  that  displacement 
of  the  stress  gage  was  identical  to  that  of  an  accderom^er  canister  at  the  same 
dq>th.  The  velocity  of  the  loading  and  unload'mg  waves  was  assumed  to  be 
constant  over  the  distance  traveled. 

In  Figure  21  the  original  positions  of  gages  at  the  upper  and  lower  depths 
are  denoted  by  HI  and  H2,  respectively.  The  displaced  positions  at  TOA  of 
the  relief  wave  are  denoted  by  and  H2(T3),  respectively.  The  location 

of  the  bottom  of  the  target  is  denoted  by  H3.  In  Figure  22  the  initial  TOA  of 
the  shock  wave  at  each  gage  d^th  is  labeled  as  77  and  72.  The  TOA  of  the 
relief  wave  at  the  lower  and  upper  depths  is  labeled  T3  and  T4,  respectively. 
Given  the  time  for  a  shock  wave  to  travel  these  two  known  distances,  and 
knowing  the  velocity  over  one  of  those  distances,  the  velocity  over  the  second 
distance  may  be  determined.  At  the  upper  dq;)th  and  lower  dq)ths  (respec- 
tively),  the  difference  in  time  between  loading  and  unloading  wave  arrival  is 
given  by: 


fTA  _  -  HI  ^  H3  -  H1(T4) 

Cl  Q 

(72  -  72)  -  ^  H3  -  772(72) 


whwe 


Q  =  the  loading  wave  velocity 
C(/  =  the  unloading  wave  velocity 

These  equations  may  be  rearranged  to  solve  for  the  unloading  wave  velocity. 
For  gages  located  at  the  upper  dq)th. 


H3  -  H1{T4) 


(T4  -  Tl)  - 


H3  -  HI 
Cl 


(1) 


For  gages  located  at  the  lower  depth,  the  unloading  wave  velocity  may  be 
determined  from 
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H3  -  H2{T3) 


{J3  -T2)  - 


H3  ~  m 


(2) 


Equation  (1)  was  applied  to  the  stress  records  at  the  upper  dq>th  and  Equa¬ 
tion  Q)  was  applied  to  the  stress  records  at  the  lower  dq)th  for  Test  22.  The 
results  are  presented  in  Table  6.  The  loading  wave  speed,  Q,  used  in  these 
calculations  was  that  determined  by  the  linear  fit  to  the  TOA  data;  i.e., 

332  mps  (see  Figure  20).  The  thickness  of  the  gage  was  also  considered  in 
these  ^culations. 


1  Table  6  I 

1  Unloading  Wave  Speed,  Cq.  as  Determined  by  the  Analytic  Method  for  1 
1  Test  22.  Assuming  a  Constant  Loading  Wave  Speed  | 

Oege  Name 

Depth,  cfii 

Cu.  "»P*  1 

HRSE-1 

15.2 

.,0 

NC8S-1 

15.2 

670  1 

CCBS-1 

15.2 

777  1 

1  HR$E-2 

30.5 

506  1 

1  NC8S-2 

30.5 

555  1 

The  variation  in  the  determined  unloading  wave  speed  may  be  attributed  to 
experiment  error  in  placing  gages  in  the  testbed,  or  Ae  assumption  that  the 
displacement  of  the  stress  gage  is  identical  to  that  of  an  adjacent  accelerometn 
canistm:.  However,  the  most  plausible  explanation  is  perhaps  the  error  in  as¬ 
suming  the  loading  wave  travels  at  constant  velocity  t^ugb  die  dqitb  of  the 
testbed.  This  may  be  seen  in  the  values  of  the  Cy  determined  by  the  stress 
gages  at  the  lower  depth.  If  a  slower  loading  wave  velocity  is  chosen  for  the 
dqtths  where  stresses  are  lower,  the  determined  unloading  wave  speed  will  be 
larger. 

This  premise  was  investigated  for  the  data  of  Test  22.  A  loading  wave 
speed,  Qj,  was  determined  for  that  portion  of  the  testbed  between  the  two 
instrumented  dqiths  by  considering  only  the  TOA  data  at  the  two  d^tbs.  A 
different  loading  wave  speed,  Q^,  was  determined  by  the  linear  fit  tfoough  the 
TOA  data  at,  and  below,  the  lower  instrumented  dqpth.  For  gages  at  the  up¬ 
per  d^th,  the  mathematical  description  of  the  time  required  for  the  shock 
wave  to  make  the  round  trip  is. 
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(TV  -  77) 


H2  -  HI  H3  ~  H2 
Cu  * 


H3  -  HHJ4) 
Cu 


where 

Qi  =  the  loading  wave  speed  between  the  upper  and  lower  instrumented 
depths 

Ci2  -  die  loading  wave  speed  between  the  lower  depth  and  the  bottom  of 
the  sample 

This  may  be  rearranged  to  solve  for 


Cv 


m  -  H1(T4) 


{T4  -  Tl)  - 

H2  -  HI  ^  H3  -  H2 

(3) 


Equation  (3)  is  used  to  determine  Cy  from  die  stress  records  for  the  three 
gages  at  the  upper  dqith.  Equation  (2)  is  still  applicable  for  the  two  stress 
gages  at  the  lower  depth  by  substituting  Q.  for  Cl,.  The  dashed  line  presented 
in  Figure  24  is  the  linear  fit  through  the  TOA  data  at  the  two  instrumented 
d^ths.  The  loading  wave  speed  determined  by  the  slope  of  this  line 
(3M  mps)  is  the  assumed  Q,  over  the  distance  from  HI  to  H2  within  the 
testbed.  The  solid  line  in  Figure  24  is  a  linear  fit  through  instruments  at  and 
below  the  lower  d^th.  The  wave  speed  determined  from  a  linear  fit  through 
diis  data  is  242  mps.  This  is  substituted  as  in  Equation  (3)  and  as  Q  in 
Equation  Q)  to  determine  the  unloading  wave  speed.  The  results  of  this 
anidysis  are  presented  in  Table  7. 

While  this  analysis  does  not  yield  a  unique  value  for  the  unloading  wave 
speed  in  the  testbed,  it  does  provide  insight  on  how  future  tests  may  be  instru¬ 
mented  in  order  to  determine  the  changing  loading  wave  speed  through  the 
d^th  of  the  testbed.  An  advantage  of  the  S&V  canister  for  this  type  of  analy¬ 
sis  is  that  the  stress  and  velocity  measurements  are  made  in  the  same  canister; 
dierefore  the  displacement  of  the  stress  gage  is  known  as  a  function  of  time. 
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Figure  24.  Determination  of  a  varying  loading  wave  speed  through  the  depth  of  the  target 


Table  7 

Unloading  Wave  Speed.  Cy.  as  Oatarmined  by  tha  Analytic 
Method  for  Tost  22,  Assunning  a  Varying  Loading  Wave  Spaed 


Oaoa  Nama 

Oapth,  cm 

C«.  mpa 

HRSE-1 

15.2 

1977 

NC8S-1 

15.2 

1253 

CCBS-1 

15.2 

1685 

HRSE-2 

30.5 

1389 

NCBS-2 

30.5 

1712 
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4  Analysis  of  Nonnormal 
Impacts 


Requirement 


As  mentioned  in  Section  3,  it  is  not  probable  to  have  a  completely  normal 
impact  b^ween  the  projectile  and  sample  when  testing  with  a  ^ft  diameter 
projectile.  A  full  understanding  of  the  interaction  between  the  projectile  and 
the  san^>le  is  required  in  order  to  use  the  large  diam^er  gas  gun  for  dynamic 
material  properties  testing. 

In  die  analysis  of  TOA  measurements  presented  earlier,  it  was  assumed 
that  the  proje(^e  inq)acted  the  testbed  in  a  perfectly  nom^  fashion.  In  order 
to  determine  the  propagation  velocity  of  the  shock  front,  the  TOA  at  eadi 
gage  location  was  plotted  against  the  rdative  distance  between  the  top  surfaces 
of  different  instruments  (for  the  first  eleven  tests).  For  Test  22,  the  TOA  was 
plotted  against  the  actual  dq)th  to  the  top  surface  of  each  instrument  in  the 
tesfoed.  The  wave  speed  in  the  sample  was  then  determined  by  die  inverse 
slope  of  the  linear  fit  through  this  data.  For  diis  type  of  analysis,  the  value  to 
be  plotted  as  "dqith"  is  die  perpendicular  distance  between  the  gage  and  the 
loading  shock  front.  For  a  perfecdy  normal  impact,  this  distance  is  sinqily 
the  vertical  distance  from  the  top  of  die  testbed  to  die  top  surface  of  die 
instrument.  For  a  nonnormal  imp^,  however,  the  gage  depth  is  a  fimction 
of  the  projectile  velocity,  the  angle  of  impact,  and  the  propagation  vdocity  of 
die  loading  wave. 


Description  of  Nonnormal  Impact 


Several  assumptions  are  required  for  the  analysis  of  nonnormal  impacts. 
These  are: 

•  The  surface  of  the  impacting  plate  is  flat. 

•  The  surface  of  the  testbed  is  flat. 
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•  Tbe  surface  of  the  testbed  is  perpendicular  to  the  cent^line  axis  of  die 
testbed. 

•  The  centerlines  of  the  projectile  impact  plate  and  the  testbed  are 
collinear. 

•  The  projectile  does  not  rotate  after  initial  contact  with  the  testbed. 

•  Tbe  loading  wave  velocity,  Q,  is  a  constant. 

In  a  nonnonnal  impact,  a  single  point  at  die  outer  edge  of  the  projectile 
impact  plate  strikes  the  top  surface  of  the  tesdied  first  (see  Figure  25).  The 
projectile  continues  to  be  driven  into  the  test  sample.  The  last  part  of  the 
impact  plate  to  come  into  contact  with  the  sample  is  the  point  on  the  face  of 
the  plate  diametrically  opposed  to  the  first  point  of  impact.  Tbe  time  delay, 
AT,  between  the  first  point  of  impact  and  the  last  point  of  impact  is  a  function 
of  the  projectile  velocity,  Vel^  the  projectile  diameter,  D,  and  the  impact 
angle,  0,  according  to  &e  relation. 


AT 


D  *  Sin(e) 
Vel, 


The  loading  wave  in  the  sample  emanates  from  the  first  point  of  impact, 
but  does  not  become  fully  developed  until  tbe  last  point  on  the  projeaile 
impacts  the  surface  of  the  testbed.  The  angle,  a,  at  which  the  normal  loading 
wave  propagates  through  the  testbed  is  a  function  of  tbe  propagation  velocity, 
Q,  in  tbe  testbed;  tbe  difference  in  time,  AT,  b^een  tbe  first  and  last  points 
of  inqiact  of  tbe  projectile;  and  the  angle  of  impact  (0)  of  the  projectile;  ac¬ 
cording  to  tbe  relation. 


a 


ran"‘ 


C^*AT 
D  *  Cos(e) 


Substituting  for  AT  from  above,  this  expression  becomes 


a  ■  Tan~' 


q  *  Taniff) 
VW 

F 


(4) 


From  this  equation,  it  is  evident  that  as  die  projectile  impact  velocity  in¬ 
creases,  tbe  angle  at  which  the  wave  propagates  through  the  sample  is 
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Figure  25.  Schematic  of  nonplanar  impact 


minimized,  for  a  givm  angle  of  impact,  6.  On  the  other  hand,  for  testbed 
maturials  widi  a  higher  loading  wave  speed,  the  effect  is  increased. 


Description  of  CLAP  Code 

A  computer  code  was  written  (in  BASIC)  to  assist  in  the  analysis  of  non¬ 
normal  impacts.  The  code  name,  CLAP,  is  an  acronym  for  loading  wave 
(CJ  angle  (A)  program  (P).  The  code  is  used  to  cali^ate  the  first  point  of 
impact  on  the  testbed  surface,  the  angle  of  impact  of  the  projectile,  and  the 
loading  wave  speed  in  the  testbed. 


Determination  of  First  Point  of  Impact  on  the  Surface  of  the  Testbed 
and  Angle  of  Impact 

For  Test  22,  three  piezoelectric  crystals  were  placed  on  the  surface  of  the 
testbed  to  determine  the  angle  of  projectile  impact.  The  differences  in  TOA 
for  the  crystals  were  calculated  using  the  first  crystal  impacted  as  the 
reference  for  the  remaining  crystals.  Assuming  that  the  crystals  lie  in  a  hori¬ 
zontal  plane  perpendicular  to  the  centerline  of  the  testbed,  differences  in  the 
TOA’s  for  the  t^ee  crystals  indicates  that  die  projectile  was  tilted  as  it 
inqiacted  the  testbed  (see  Figure  26).  At  the  time  the  last  crystal  was 
impacted,  the  product  of  the  impact  velocity  and  the  difference  in  TOA  yields 
the  relative  height  of  the  projectUe  (at  the  crystal  locations)  with  respect  to  the 
testbed  surface,  assuming  the  projectile  does  not  rotate  after  impact.  With  the 
centerlines  of  the  projectile  and  testbed  aligned,  die  (x,y)  coordinates  of  the 
crystals  may  be  transferred  to  the  face  of  the  projectile.  Three  points  on  the 
projectile’s  impacting  surface  are  now  defined,  with  the  origin  of  the  coordi¬ 
nate  syston  at  the  top  center  of  the  testbed  surface.  The  equation  of  the  plane 
throu^  these  three  points  is  then  determined. 

Using  the  equation  describing  the  plane  of  the  impact  plate,  the  first  point 
of  contact  onto  die  testbed  surftice  is  determined.  The  impact  plate  is  divided 
into  1440  equal  segments  ('A  degree  increments).  For  each  increment,  the  z 
coordinate  at  the  outer  radius  of  the  plate  is  calculated.  The  increment  which 
results  in  the  minimum  z  rqiresents  Ae  first  point  of  impact  at  the  surface  of 
the  testbed.  The  angle  formed  by  the  horizontal  testbed  surface  and  the  line 
from  the  first  point  of  impact  through  the  center  of  the  impact  plate  is  the 
angle  of  impact,  B  (see  Figure  26). 
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26.  Schematic  of  parameters  used  to  determine  the  first  point  of  impact  on  the  testbed  surface  for  nonplanar  impacts 


Determination  of  the  Loading  Wave  Speed  in  the  Testbed  Using 
TOA  Measurements  at  the  Surface  of  the  Testbed 

To  b^in,  the  loading  wave  speed  is  determined  by  the  procedure  described 
in  Section  3.  That  is,  a  normal  impact  is  assumed  and  the  loading  wave  speed 
is  the  inv^e  slope  of  a  linear  fit  t^ugh  the  TOA  versus  d^th  plot.  A 
correlation  coefficient  is  determined  that  r^resents  the  quality  of  the  linear  fit 
through  the  expwiment  data.  For  a  perfect  fit  this  coefficient  is  1 .  For  any 
case  odier  than  a  perfect  fit  this  coefficient  is  less  than  1.  Using  Equation  (4), 
the  angle  (a)  at  which  die  normal  loading  wave  propagates  through  the  testbed 
is  calculated.  The  'uepdi"  of  each  gage  location  is  recalculated  as  the  perpen¬ 
dicular  distance  between  the  gage  location  and  the  new,  fully-developed 
loading  wave  (with  angle  a).  This  is  illustrated  in  Figure  27. 

A  new  loading  wave  speed  is  determined  from  the  linear  fu  through  a  plot 
of  TOA  versus  (recalculated)  depth.  The  correlation  coefficient  is  also  calcu¬ 
lated.  As  the  values  of  dq>di  are  adjusted  for  each  gage  location,  the  data 
falls  n  I'e  closely  on  a  line,  and  thus  increases  the  correlation  coefficient. 
Using  the  new  value  for  the  loading  wave  velocity.  Equation  (4)  is  used  once 
again  to  determine  the  angle  at  which  the  plane  wave  propagates  through  the 
testbed.  The  ”dq)th"  of  each  gage  is  recalculated,  and  a  new  loading  wave 
velocity  is  determined  as  before.  This  itoative  process  continues  until  the 
correlation  coefficient  converges  to  a  constant  value. 

Test  22  was  the  only  test  that  used  TOA  crystals  at  the  surface  of  the 
testbed  to  determine  the  angle  of  impact.  CLAP  was  used  to  evaluate  the 
angle  of  inq)act  (0),  the  furst  point  of  impact,  the  angle  at  which  the  loading 
wave  travelled  through  the  testbed  (a),  and  the  loading  wave  velocity  {Q),  as 
described  above.  A  series  of  screen  dumps  from  the  code  are  included  as 
Figures  28  through  35.  Figure  28  is  an  introductory  screen  allowing  the  user 
to  sdect  one  of  two  methods  of  analysis.  The  first  selection  is  described 
above  and  the  second  selection  is  a  trial  and  error  method,  described  below. 
Data  required  by  the  program  (e.g.,  gage  locations,  TOAs,  ^c.)  may  be 
entwed  through  the  keyboard  or  by  a  data  file.  Figure  29  lists  the  data  used 
to  calculate  the  equation  of  the  plane  through  the  impact  plate.  Also  presented 
are  the  coefficients  for  the  equation  of  the  plane.  Figure  30  indicates  the  first 
point  of  impact  of  die  projec^e  onto  the  testbed  surface.  Figure  31  lists  the 
coordinates  of  the  gage  locations  in  the  testbed.  This  coordinate  system 
defines  the  origin  at  the  top  cento'  of  the  testbed,  with  the  positive  x-axis 
pointing  due  east  when  looking  at  a  plan  view  of  the  testbed  (see  Figure  3). 
The  values  for  z  indicate  the  d^th  to  the  top  surface  of  the  gage  (note:  this  is 
a  left  handed  coordinate  system).  Figure  32  lists  the  least  squares  curve  fit 
data  for  the  data  set  of  Figure  31 .  The  plot  of  Figure  33  is  identical  to  the 
TOA  versus  d^th  plot  of  Figure  20.  Figure  34  lists  the  least  squares  curve 
fit  data  after  one  iteration  on  the  angle  at  whidi  the  loading  wave  propagates. 
The  plot  of  Figure  35  shows  the  TOA  data  plotted  against  the  recalculated 
values  for  the  dqith  of  the  gages.  The  loading  wave  speed  changed  from 
332  mps  to  336  mps  after  one  iteration,  and  the  correlation  coefficient 
decreased  from  0.98222  to  0.97923.  Further  iterations  do  not  significantly 


Chapter  4  Analyaia  of  Nonplanar  Impacts 


CENTER  OF 
IMPACT  PLATE 


44 


Chapter  4  Analysis  of  Nonplanar  Impacts 


Figure  27.  Orientation  of  coordinate  system  used  to  determine  the  loading  wave  speed 


TOA  CRYSTAL 
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Figure  28.  Elevation  view  indicating  "depths”  used  in  determining  the  loading  wave  speed 


Program  CLAP  . . .  Loading  Wave  (Cl)  velocity  and  Angle  Program 
February  09,  1993,  program  version  2.3 


US  Army  Corp  of  Engineers 
Waterways  Experiment  Station 
Structures  Laboratory,  Explosion  Effects  Division 
Vicksburg,  Mississippi 


1  -•  Calculate  Cl  by  a  linear  regression  curve  fit  to  data 

2  --  Calculate  Cl  by  fitting  a  plane  to  the  wave  front 


Enter  '1',  '2'  or  'q'  to  quit 


Figure  29.  Introductory  screen  of  the  CLAP  code 
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I,  .  I-  I—  ■  Entered  Impact  Plate  Data  -  "■  . . — 

The  following  data  has  been  entered 
Test  name:  test  22 

Location  of  TOA  gage  1  -  (  0  ,  .1524  ) 

Location  of  TOA  gage  2  -  (  .32385  ,-.32385  ) 

Location  of  TOA  gage  3  -  (-.32385  ,-.32385  ) 

Time  of  arrival  for  TOA  gage  1  -  558.247  msec 
Time  of  arrival  for  TOA  gage  2  -  558.246  msec 
Time  of  arrival  for  TOA  gage  3  -  558.145  msec 
Projectile  velocity  -  50  m/sec 

The  general  equation  for  the  equation  of  a  plane  is  Ax+By+Cz+D-0 

The  following  coefficients  have  been  calculated  for  the  projectile 
A  — 2.402651E-03 
B  — 1.668394E-03 
C  -  .3084671 
D  —8.609789 


strike  'c'  to  continue  or  'q'  to  quit 


Figure  30.  CLAP  screen  listing  data  used  for  determining  the  equation  of  the  plane  of  the 
projectile  and  the  calculated  coefficients  of  the  plane 
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li—i— . .  "■  First  Point  of  Impact  »  ' 

Using  the  equation  for  the  plane  representing  the  impact  plate 
previously  calculated,  the  first  point  of  impact  on  this  plate 
will  now  be  determined. 

Test  name:  test  22 

Entered  Impact  plate  diameter  -  47  inches 

Entered  impact  plate  diameter  -  1.1938  meters 

Entered  projectile  velocity  ••  50  m/sec 


The  first  point  of  impact  occurs  on  the  impact  plate  at  236.25  degrees. 
This  occurs  at  the  plate  radius  of  0.5969  meters. 

The  maximtun  angle  of  Impact  is  0.543  degrees.  (  9.48  mrad  ) 


Strike  'c'  to  continue  or  'q'  to  quit 


Figure  31.  CLAP  screen  indicating  the  first  point  of  impact  onto  the  surface  of  the  testbed  and 
the  angie  of  impact 
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Figure  32.  CLAP  screen  listing  the  data  used  to  determine  the  initial  loading  wave  speed  (with 
no  iterations  in  CLAP) 
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Least  Square  Curve  Fit 


Test  ID  :  Test  22,  All  gages  Included 


NuDd>er  of  Iterations 
Wave  Propagation  Angle 
Slope 

Cl  based  on  slope 
Intercept 
Standard  Error 
Correlation  Coefficient  (R) 


0  degrees 

3.011879E-03 

332.0187  meter/sec 

.558178 

7.85667E-05 

.9822238 


Coefficient  Determination  (R2)  ••  .9647636 


•i—  strike  '1'  to  iterate,  'p'  to  plot,  or  'q'  to  quit  or  'c'  to  continue 
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Figure  34.  CLAP  screen  showing  TOA  vs.  depth  data  (prior  to  iterating  with  CLAP> 


Least  Square  Curve  Pit 


Test  ID  ;  Test  22,  All  gages  Included 


Number  of  iterations  - 
Wave  Propagation  Angle  •• 
Slope  - 
Cl  based  on  slope  ~ 
Intercept  ~ 
Standard  Error  - 
Correlation  Coefficient  (R)  - 
Coefficient  Determination  (R2)  <- 


3.602976  degrees 

2.993212E-03 

334.0893  meter/sec 

.558279 

8.979428E-05 

.9767154 

.9539731 


—  strike  '1'  to  iterate,  'p'  to  plot,  or  'q'  to  quit  or  'c'  to  continue 
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Figure  36.  CLAP  screen  showing  TOA  vs.  depth  data  (after  a  single  iteration  with  CLAP) 


change  the  loading  wave  speed  or  the  correlation  coefficient.  For  example, 
after  an  additional  four  iterations  the  d^ennined  loading  wave  speed  is 
336  mps  and  die  correlation  coefficient  is  0.97913.  The  angle  at  which  the 
wave  propagated  through  the  testbed  was  dttermined  to  be  3.6S  degrees. 

As  a  result  of  this  analysis,  the  stated  locations  of  several  gages  in  die 
testbed  appear  to  contain  errors.  This  is  because  the  correlation  coefficient 
should  increase  after  each  iteration,  and  for  this  data  set  the  coefficient  de¬ 
creased.  Several  factors  could  result  in  inaccurate  gage  locations:  1)  gages 
were  not  placed  properly  in  the  testbed,  2)  gages  locations  were  not  measured 
or  recorded  properly,  or  3)  the  gage  locations  changed  when  the  testbed  was 
relocated  from  foe  prqiaration  area  to  foe  gas  gun  facility.  An  additional 
source  of  error  could  simply  be  foe  tolerance  at  which  foe  gages  were  placed 
in  foe  testbed  (±3  mm). 

From  foe  plot  of  Figure  35,  it  is  noted  that,  with  foe  exception  of  five 
TOA  measurements,  foe  data  set  is  quite  linear.  These  five  measurements  are 
indicated  as  "oudiers"  in  foe  figure.  As  an  exercise,  CLAP  was  used  to 
evaluate  a  modified  data  set  that  removed  these  measurements.  A  plot  of  this 
data  set  after  a  single  iteration  is  shown  in  Figure  36.  The  loading  wave 
speed  determined  using  this  reduced  data  set  was  340  mps.  The  angle  at 
which  foe  wave  propagated  through  foe  sample  was  determined  to  be 
3.71  degrees,  llie  correlation  coefficient  increased  to  0.99766,  indicating  a 
tighter  data  set  than  that  of  Figure  35. 

By  accounting  for  nonnormal  impact,  foe  loading  wave  speed  increased  by 
0.6  percent  when  considering  foe  entire  data  set,  and  increased  by  3  percent 
when  considering  a  modified  data  set.  This  difference  is  not  significant  in 
either  case.  What  is  significant,  however,  is  that  this  analysis  provides  an 
ind^endent  check  on  foe  location  of  foe  gages.  It  will  also  assist  in  determin¬ 
ing  foe  required  placement  tolerances  on  future  tests. 


Determination  of  the  Loading  Wave  Speed  in  the  Testbed 
on  a  Trial  and  Error  Basis 

If  foe  angle  of  impact  of  foe  projectile  is  not  known  or  cannot  be  d^er- 
mined,  foe  CLAP  code  uses  a  trial  and  error  method  for  determining  foe 
loading  wave  speed  in  foe  soil.  This  method  tries  different  combinations  of 
first  points  of  impact  and  loading  wave  propagation  angles  (a)  to  determine 
which  condition  will  produce  foe  highest  correlation  coefficient  from  foe  linear 
regression  curve  fit. 

To  do  this,  CLAP  sets  foe  first  point  of  impact  at  0  degrees.  Then  foe 
angle  between  foe  loading  wave  and  foe  testbed  surface  is  varied  from  0  to 
12  degrees  in  14-degree  increments.  For  each  increment,  foe  perpendicular 
distance  (or  depth)  between  foe  gage  location  and  foe  loading  wave  plane  is 
calculated,  as  described  above.  A  linear  fit  through  foe  TOA-versus-depth 
data  is  made  and  foe  correlation  coefficient  is  recorded.  The  first  point  of 
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Figure  37.  CLAP  screen  showing  TOA  vs.  depth  data  for  a  modified  input  data  set  (after  a  single  iteration  with  CLAP) 


impact  is  dien  incremoited  by  1  d^ee.  As  before,  the  angle  between  the 
loading  wave  and  the  testbed  surface  in  incremented.  A  linear  fit  is  made  and 
die  corrdation  coefficient  for  each  increment  is  calculated.  This  procedure  is 
continued  as  the  first  point  of  impact  varies  from  0  to  3S9.7S  degrees.  The 
combination  of  the  first  point  of  impact  and  the  angle  of  the  loading  wave 
which  produced  die  best  (or  highest)  correlation  coefficient  is  said  to  r^resent 
die  actual  test  condition.  The  invose  slope  of  the  linear  fit  dirough  the  data 
for  diis  combination  is  reported  as  the  loading  wave  speed. 

The  graphical  procedure  for  determining  the  loading  wave  speed,  pressed 
in  Section  3,  is  very  sensitive  to  the  accuracy  of  the  gage  placement.  The 
differrace  in  hei^t  between  the  top  surface  of  some  gage  types  (at  the  same 
level)  is  on  die  same  order  as  the  placement  tolerances.  A  comparison  of 
values  for  the  loading  wave  speed  as  dttermined  by  the  grsqihic^  procedure, 
vmsus  those  determined  by  the  CLAP  trial  and  error  procedure,  can  better 
describe  actual  test  conditions. 

Agreement  b^een  the  loading  wave  speed  determined  by  the  CLAP  code 
and  the  value  determined  by  the  gri^hical  procedure  (discussed  in  Section  3) 
would  indicate  normal  impact.  However,  large  values  for  alpha  would  not 
necessarily  indicate  a  large  impact  angle,  since  the  angle  at  which  the  wave 
propagates  through  the  sample  is  a  function  of  the  projectile  velocity. 

Repeated  in  Table  8  is  the  value  of  the  loading  wave  speed  determined  by 
the  graphical  method  presented  in  Table  5,  Seaion  3.  The  trial  and  error 
procedure  within  CLAP  was  used  to  d&ermine  the  loading  wave  speed  and 
angle  of  propagation  for  the  tests  where  foe  angle  of  impact  could  not  be 
determined  (all  ex^t  Test  22).  The  results  are  presented  in  Table  8.  The 
correlation  coefficient  for  each  method  is  listed  in  foe  table.  Note  that  in  each 
instance  the  wave  speed  decreases  or  remains  constant  and  foe  correlation 
coefficient  increases.  Using  Equation  (4),  Section  4,  foe  angle  of  impact  (9) 
was  calculated  to  check  foe  consistency  of  foe  values  for  Q  and  a  determined 
using  CLAP.  The  values  presented  in  Table  8  are  consistent  and  similar  to 
foe  value  measured  on  Test  22.  Hence  foe  trial  and  error  method  used  within 
CLAP  provides  a  mechanism  to  analyze  foe  geometry  of  impact  without  sur¬ 
face  TOA  measurements. 
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Tables 

Comparison  of  Loading  Wave  Speed  Determined  by  Graphical  Pro¬ 
cedure  and  by  Trial  and  Error  CLAP  Calculation 


TMt 

No. 

Sand 

Typa 

QrapMeal  Maihod 

CLAI^,  Trial  arwi  Error  Method  | 

c,. 

mpa 

Corralaiion 

Coaffidant 

Ci. 

mpa 

— 

Ccrrslation 

Coaffidant 

§. 

mlliradiana 

e 

MMonry 

399 

0.93697 

363 

0.95152 

3.99 

5.4 

7 

Mmoiwv 

491 

0.98723 

465 

0.99941 

3.24 

4.5 

8 

Masonry 

S03 

0.85381 

506 

0.93056 

8.65 

10.7  1 

10 

Masonry 

328 

0.98176 

275 

0.99986 

4.73 

14.2  1 

11 

Masonry 

417 

0.96161 

341 

1.0 

5.72 

14.8  1 

12 

Masonry 

398 

0.98081 

349 

1.0 

10.09 

21.3  1 

13 

Masonry 

337 

0.95557 

333 

0.99608 

5.72 

15.9  1 

1  14 

Masonry 

323 

0.94881 

296 

0.96274 

3.74 

12.3 

1  15 

Masonry 

323 

0.89308 

288 

0.98996 

7.44 

26.7 

1  15 

Socorro 

280 

0.95149 

276 

0.97939 

8.65 

25.5 

1  1"^ 

Socorro 

309 

0.97052 

270 

0.97786 

1.50 

5-3  1 
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5  Conclusions  and 
Recommendations 


Conclusions 

Twelve  tests  were  conducted  with  the  DNA/WES  Ground  Motion  Test 
Facility  with  sand  testbeds  as  targets.  Nine  tests  used  Masonry  sand  and  three 
tests  used  Socorro  Plaster  sand  as  the  target  testbed  material.  A  variety  of 
stress  and  motion  transducers  were  included  in  each  testbed  to  investigate  the 
environment  in  the  testbeds,  and  to  determine  the  robustness  of  the  various 
instruments. 

These  experiments  provided  valuable  information  on  the  mechanics  of 
testing  with  the  gim.  Analysis  of  the  data  provided  insight  for  instrumentation 
requirements  and  placement  procedures/tolerances  for  future  tests.  The  results 
indicate  that  the  facility  is  a  viable  tool  for  conducting  dynamic  tests  for 
instrumentation  development  and  material  properties  studies. 


Recommendations 

Future  tests  with  the  gun  will  incorporate  more  diagnostic  measurements  in 
the  testbed.  These  will  include  time-of-arrival  (TO A)  measurements  at  the 
surface  of  the  testbed  to  determine  the  velocity  and  angle  of  impact  of  the 
projectile.  The  TOA  of  the  shock  wave  at  various  depths  in  the  sample  will 
be  measured  (as  was  done  for  Test  22)  to  supplement  the  TOA  measured  at 
the  active  gage  locations.  The  inclusion  of  additional  TOA  measurements  will 
allow  for  more  exact  determination  of  the  propagation  velocity  of  the  loading 
and  unloading  waves  in  the  sample.  Also,  the  S&V  canister  will  be  fielded  to 
allow  for  more  precise  determination  of  the  displacement  of  a  stress  sensor, 
which  will  in  turn  allow  a  more  definitive  calculation  of  the  unloading  wave 
speed. 

In  order  to  make  optimum  use  of  the  gun  to  validate  or  test  ground  motion 
sensors  in  granular  materials,  a  testbed  material  is  required  that  will  have 
consistent  and  predictable  properties.  WES  is  investigating  the  selection  of  a 
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standard  granular  material  (perhaps  synthetic)  and  the  development  of 
fnnplarj»jw<»nt  procedures  which  will  provide  a  testbed  with  well-controlled  and 

consistent  properties. 
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Refarancas 


Appendix  A 
Notation 


Q  Average  loading  wave  velocity  in  a  testbed  (velocity  of  distur¬ 
bance  caused  by  projectile  impact) 

Cu  Loading  wave  velocity  between  the  two  instrumented  depths  in  a 
testbed 

Qj  Loading  wave  velocity  between  the  lower  instrumented  d^th  and 
bottom  of  a  testbed 

Cy  Unloading  wave  velocity  in  a  testbed  (wave  resulting  from  reflec¬ 
tion  from  bottom  of  tested) 

D  Projectile  diameter 

HI  Upper  instrumented  depth  in  a  testbed 

H1(T3)  Location  of  an  instrument  originally  in  the  upper  depth,  at  TOA 
of  relief  wave 

H2  Lower  instrumented  d^th  in  a  testbed 

H2(T4)  Location  of  an  instrument  originally  in  the  lower  depth,  at  TOA 
of  relief  wave 

TOA  Time-of-arrival  (of  shock  waves) 

77  Initial  TOA  at  a  gage  in  the  upper  depth 

72  Initial  TOA  at  a  gage  in  the  lower  d^th 

T3  TOA  of  the  relief  wave  at  a  gage  in  the  lower  dqjth 

T4  TOA  of  the  relief  wave  at  a  gage  in  the  upper  d^th 
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AT  Difference  in  time  between  the  first  and  last  points  of  impact  of 
the  projectile  onto  the  surface  of  a  testbed 

Velf  Projectile  velocity 

ct  angle  at  which  the  loading  wave  propagates  through  the  testbed 

9  angle  at  which  the  projectile  impacts  the  testbed 
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Table  B1  (Continued) 


Depth  B  16.2  cm 
lUdkM  •  30.6  cm 


Depth  •  30.6  CM 
Retiiue  «  16.2  CM 


(Shmmt  3  of  4! 


'  Located  in  center  of  testbed  (radius  ^  0  cm) 
’  Oriented  horizontally 
’  No  data 
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Figure  C2.  Wet  density  and  moisture  content  profiles  for  Test  7 
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Figure  C4.  Wet  density  and  moisture  content  profiles  for  Test  10 
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Figure  C7.  Wet  density  and  moisture  content  profiies  for  Test  13 
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Figure  C8.  Wet  density  and  moisture  content  profiles  for  Test  14 
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Figure  C9.  Wet  density  and  moisture  content  profiles  for  Test  15 


Cl 


A 


Appendix  C  Wet  Density  end  Moisture  Content  Profiles  for  Tests  in  Send 


Figure  C10.  Wet  density  and  moisture  content  profiies  for  Test  16 
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Fifsure  Cl  1 .  Wet  density  and  moisture  content  profiies  for  Test  1 7 
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Figure  Cl  2.  Wet  density  and  moisture  content  profiles  for  Test  22 
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Figure  D1 .  Stress  and  impulse  wave  forms  from  Test  6  (Continued) 
D2 
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IMPULSE  MPAXS 


Figure  D2.  (Concluded) 
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Figure  D5.  Stress  and  impulse  wave  forms  from  Test  1 1  (Continued) 
DIO 
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Figure  D5.  (Concluded) 
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Figure  D6.  Stress  and  impulse  wave  forms  from  Test  12  (Continued) 


Appendix  D  Active  Measurement  Data  for  Tasts  in  Sand 


Figure  D7.  Stress  and  impulse  wave  forms  from  Test  13  (Continued) 
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Figure  D1 1 .  (Concluded) 
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Figure  D12.  Stress  and  impulse  wave  forms  from  Test  22  (Continued) 
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Figure  D13.  Acceleration,  velocity,  and  displacement  wave  forms  from  Test  6  (Sheet  1  of  3) 
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Figure  D14.  Acceleration,  veiocity,  and  displacement  wave  forms  from  Test  7  (Sheet  1  of  3) 
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Figure  D14.  (Sheet  2  of  3) 
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Figure  D15.  Acceleration,  velocity,  and  displacement  wave  forms  from  Test  8  (Sheet  1  of  3) 
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Figure  D15.  (Sheet  3  of  3) 
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Figure  D16.  Acceleration,  velocity,  and  displacement  wave  forms  from  Test  10  (Sheet  1  of  3) 
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Figure  D16.  (Sheet  3  of  3) 
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Figure  Dl7.  Acceleration,  velocity,  and  displacement  wave  forms  from  Test  11  (Sheet  1  of  3) 
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Figure  D17.  (Sheet  3  of  3) 
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Figure  D18.  Acceleration,  velocity,  and  displacement  wave  forms  from  Test  12  (Sheet  1  of  3) 
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Figure  D18.  (Sheet  3  of  3) 


Appendix  D 


Active  Measurement  Data  for  Tests  in  Sand 


ACCELERATION  -  G-S 


658 . 6  659 . 1  659 . 6  660 . 1  660 . 6  661 . 1  661.6  662 . 1  662 . 6  663 . 1  663 . 6 

TIME  -  MSEC 


Figure  D19.  Acceleration,  velocity,  and  dispiacement  wave  forms  from  Test  13  (Sheet  1  of  3) 
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Figure  D20.  Acceleration,  velocity,  and  displacement  wave  forms  from  Test  14  (Sheet  1  of  3) 
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Figure  D22.  Acceleration,  velocity,  and  displacement  wave  forms  from  Test  16  (Sheet  1  of  3) 
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Figure  D22.  (Sheet  3  of  3) 
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Figure  D23.  (Sheet  4  of  4) 
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Figure  D24.  Acceleration,  velocity,  and  displacement  wave  forms  from  Test  22  (Sheet  1  of  3) 
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